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We show that a spatially-dependent magnetic field can induce a persistent 
spin current in a spin-1/2 Heisenberg antiferromagnetic ring, proportional to 
the solid angle subtended by the magnetic field on a unit sphere. The result 
is a direct consequence of Berry "parallel transport" in space. The magnitude 
of the spin current is determined by the ratio of longitudinal and transverse 
exchange interactions J\\/J± and by the magnetic field. For large magnetic 
fields the Zeeman energy strongly renormalizes the Ising term giving rise to 
a maximum spin current. In the limit of J\\/J± <C 1 the amplitude of the 
current behaves like 1 — (J\\/J±) 2 - In the opposite limit ttJ± > Jy > J± the 
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amplitude scales as J J\\/ J± exp (—Jj\\/J±)- 
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The generation of persistent currents in a ring is a consequence of the sensitivity of 
a wavefunction to its boundary conditions. The basic idea of persistent currents was put 
forward a few decades ago [1-3]. The effect of a magnetic flux threading a metallic ring, such 
that the magnetic field at the ring vanishes, can be replaced by a twist of the wavefunction 
boundary conditions. The change in these boundary conditions gives rise to an asymmetry 
between the number of electrons with K > Kp and K < —Kp, where Kp denotes Fermi 
momentum. As a result, an electric current is induced in mesoscopic ring structures. For 
large ring structures, in the thermodynamic limit, the current vanishes; in agreement with 
the fact that no spontaneous broken chirality occurs in non-interacting systems. A new 
interest in this problem has arisen for non-Fermi liquids. In particular, for Luttinger liquids 
novel behavior has recently been proposed [4-7]. 

In the last ten years the "parallel transport" idea introduced by Berry [8] for adiabatic 
processes has provided new insight into currents induced by topological effects. A typical 
situation considered by Berry was a spin-1/2 electron in a time periodic magnetic field. 
Assuming adiabatic motion Berry showed that when an electron returns to its initial position 
the wavefunction accumulates a nontrivial phase equal to the angle of rotation of the spin-1/2 
angular momentum on a sphere. This is the Berry phase [8-10]. 

Recently it has been shown [11] that for an integer spin S = 1 ferromagnet the Berry 
phase combined with the spin-orbit effect can give rise to a spin current and a corresponding 
electric field. Here we propose a new way of generating a spin current, when a spatially 
periodic magnetic field acts on a spin-1/2 antiferromagnetic ring a Berry phase in space is 
induced giving rise to a twist on the wavefunction boundary condition. As a consequence a 
persistent spin current is induced. This result can be understood within the Berry parallel 
transport idea since a spin-1/2 antiferromagnet is Lorentz invariant. 

We consider a spin-1/2 antiferromagnetic Heisenberg ring of circumference L in the 

— # 

presence of a spatially-dependent magnetic field of the "crown shape" form, h(x) = 
ho(sm0(x) cosc/)(x), sm9(x) sm<f)(x), cos9(x)), where 9{x) and (f>(x) are polar — n < <p(x) < 
7r, and azimuthal < 9(x) < 9 < n/2 angles such that <p(x) = 4>(x + L), see Fig. 1. We 
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restrict 9{x) to the "northern" hemisphere which obeys 9(x) — 9{x + L). The Hamiltonian 
for the spin-1/2 antiferromagnet in an external magnetic field is: 

" X X X 

where J± and J\\ are the transverse and longitudinal exchange coupling constants. 

We rotate S(x) in the direction of the magnetic field h(x) using the following unitary- 
transformation: 

U = exp (-^3) exp (-^2) , (2) 

where the circumference of the ring is L = Na and a 2 and a 3 are Pauli spin-1/2 matrices. 
As a result of the unitary transformation Eq. (2), the Hamiltonian H is transformed to 
H = U^HU . Due to the periodicity of the angles 9{x) and <p(x) the transformed wavefunction 
obeys the same boundary condition as the original wavefunction. 
The transformed Hamiltonian H takes the form, 

j N-l 

H-'^-T. (S+(x)e iA{x > x+1) S-(x + 1) + S-(x)e- tA ( x ' x+1) S+(x + 1)) 



N-l 

+J|, E S 3 (x)S 3 (x + 1) + 9l i B ho E S 3 (x). (3) 

x=0 x 

In Eq. (3) x + 1) is the Berry connection: 

A(x,x + 1) = = ^(1 + cos6(x))j?-dx. (4) 

To obtain Eq. (4) we have restricted 9{x) to the range < 6(x) < 0q < n/ '4, where ^0 
denotes the maximum tilt angle of the "crown-shaped" magnetic field (Fig. 1). The Berry 
connection results from rotating the coordinate of the spin-1/2 operator. Moving around 
the ring, a phase proportional to the subtended solid angle accumulates. 

Due to the unitary transformation the exchange couplings J±, J\\ are replaced by 
J± and J||, where J\\ = J\\ (cos 2 9(x)), J± = J±((l + (cos 2 9(x)))/2) with (cos 2 9(x)} = 
(1/9 q )JS° cos 2 9d9 = (1/2)(1 + (sin 29 /29 )). 
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Following Ref. [12] we use a Jordan- Wigner transformation to replace the spin-1/2 S(x) 
operator by the spinless fermions c\ c(x): 

S_(x) = e^^o^cix), S+(x) = [S-(x)}\ S 3 (x) = c\x)c(x)-\/2, n(x) = c\x)c(x). 

(5) 

It is convenient to replace the fermion c(x) by d(x), c(x) = e lF ^d(x), c\x) = d) ' (x)e~ %v<yX \ 
In addition, we choose A(x, x+l)+n(n(x) — l) = T(x) — T(x+1). As a result the Hamiltonian 
in Eq. (3) is replaced by: 

N-l v- 1 

I //(.c ) - — I I //(.c 4- 1 i - 



t N-l N-l , y\ / 

H = - -± £ (d + (x)d(x + 1) + H.c.) + J,, £ Us) - - ) (n(s + 1) 



+#/^o £ ~ \) ■ ( 6 ) 

To obtain the form given in Eq. (6), we require that the phase T(x) obeys the condition 

e i[ J 4(x,x+l)+7r(ra(x)-l)-r(x)+r(x+l)] _ j ^) 

The periodicity of the spin operator S(x) = S(x + L) together with the condition on the 
Berry connection gives rise to a twisted boundary condition for the d(x), cft(x) fermions: 

d( X + N) = -e^^=o A(x+j0+J+l)j2nN Fe -in(N+l) d ( x ^ ^ 

where Np = J2x=o n ( x ) an d L = Na. Replacing the sum in Eq. (8a) by a continuous 
integral we find, using Stokes theorem, that 

J2A(x + J,x + l + J)~f = J o J f— ^Jd^ = 7r(l-cos0o) = n(0 o ). (86) 



j=o 



Here, fl(9 ) measures the twist on the boundary conditions. Consequently, the equation for 
d(x) is replaced by 

d(x + L) = -e i * (n(9o)/,r - 1) d(a;)e i ' rJV = -e^ 1 ^), (8c) 

where 5 = Q(9o)/n + Ajv, = if N is even and An — 1 if N is odd. Equation (8c) 
expresses the twist of the boundary condition. Consequently, the "twist" of the boundary 



condition is equal to the solid angle £l(6 ), see Fig. 1. Equations (8b) and (8c) represent 
one of our major results: a geometrical rotation of the spin-1/2 coordinate gives rise to a 
Berry phase which changes the boundary conditions for the wavefunction. 

We construct the continuum representation of the Hamiltonian H. The fermion operator 
d(x) is replaced by ip(x): 

d(x) = yfal){x) = ^[e lk >e- l ^ /L) ^ 1)x R{x) + e~ ik >e^ 2w/L){s ' 1)x L(x), (9) 

where k F = -w/2a and R(x) and L(x) are the right and left moving operators at half filling: 

R( x j = 1 c ia Rc i(2ir/L)(N R -l/2)x c iVA^e R ^ 

\f2~na 

L( X ) = -^ e -^L e -i(2n/L)(N L -l/2) Xe -iV4^h_ ( 1Q ) 

y/2ira 

The representation given in Eq. (10) has been introduced in Refs. [7,13]. R and Ol 
are the particle hole excitations (non-zero modes), N R and Nl are the zero mode charge 
operators. Nr and Nl measure the additional charge with respect to the half-filled spinless 
Fermi sea. cvr and represent the zero mode coordinates which are canonical conjugates 
to Nr and Nl, respectively: 

[<XR,N R ] = [-a L ,N L ]=i, [a,J} = 2i, \J3, Q] = 2i, J = N L - N R , Q = N L + N R , 

a = a R + a L , (3 = a L - a R . (11) 
The complete set of zero mode eigenfunctions are given by 

(a R \N R ) = -}=e ia « N «, (a L \N L ) = (12a) 

V 2it v 2tt 

or 

(a\J) = ^=e iaJI \ (P\Q) = -^eW 2 , (126) 



where N R \N R ) = N R \N R ), N L \N L ) = N L N L ), J\J) = J\J) and Q\Q) = Q\Q). 

We perform a unitary transformation which removes the chiral phase (2n/L)(5 — l)x from 
Eq. (9). The Fermi momentum k° F = n/2a is shifted to k F ^ k° F such that the vacuum is 



shifted as a function of the magnetic field ho- We introduce the new unitary transformation 

V = e -(i/ 2 W-l)-q]aR e -(i/2)[(S-l)+q]a L _ ^3) 

As a result of the unitary transformation i]){x) is replaced by 

Vty(x)V = e i{k r + * q/L)x R{x) + e- i{ - k F + * q IV x L(x). (14) 

The value of q can be chosen such that the vacuum is shifted and the term gusho XXn(x) — 
1/2) is canceled [14]. 

The effect of the unitary transformation given in Eq. (13) is to change the boundary 
condition for the eigenfunction in Eqs. (12), 

V^\J)\Q) = \J)\Q), (15a) 

such that (a + An\ J) = (a\J) is replaced by 

(a + Air\J) = e i{& - 1)2 *(a\J). (156) 

Similarly we have for (f3\Q), 

(f3 + 4n\Q) = e^-^([3\Q), (15c) 

where [q] is the closest integer to q. 

Due to the boundary condition in Eq. (15b), the eigenvalues J are given by J = n/2 + 
(5 - l)/2, n = 0, ±1, ±2, ... and Q = m/2 + {q-[q])/2, m = 0, ±1, ±2, .... Using the unitary 
transformation given in Eq. (13) and the Bosonic representation in Eq. (10), we obtain the 
Bosonic form of the Hamiltonian: 

H = V^HV = V^U^HUV, H = H n=0 + H n ^ + AE , (16a) 

where H n= o is the zero mode Hamiltonian 

H n=0 = T (N 2 R + N 2 L + 2-tNrNl), (166) 
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where v = J±a(l + J\\/irJ±) and 7 = 2J\\a/nv. H n ^o is the non-zero mode Bosonic Hamil- 
tonian [14] 

H n ^ = J dxl- K(d x <pf + —(d x 6f + A cos(2a + y/TE^O + —Qx + A(K F - K° F )x) , 

(16c) 

where v = v^{\ - 7 2 ), X = ^/(l - 7 )/(l + 7), A = A/a 2 , * = ^!|a/27r 2 with = # L - fl , 
6 = 6l + Or. Here g is chosen such that the linear term gfisho is canceled and the ground 
state energy is shifted by, 



2a J ± l + 5Jj|/7rJ ± ' v * ^ L J ± (l + 5 J||/ttJ ± ) ' 

(17) 

Equation (17) shows that the shift in the Fermi momentum depends on the ratio between 
the Zeeman energy to the spin liquid lowest state excitation energy. 

The spin current is obtained from the Hamiltonian given in Eq. (16). The Heisenberg 
equation of motion for the zero modes introduced in Ref. [13] allows us to construct the 
spin current operator. We replace the electric charge by the magnetic charge s = g\in- As 
a result the spin current is given by the time derivative of the zero mode coordinate: 

s da s r - sv S/ * »-\ \ J±a ( J\\\ . . 

h - = ^ = w^ H ^-^- N ^ "^jri'-J-J- (18) 

From Eq. (18) we see that the spin current is determined by the spin liquid velocity v s = 
J±a/h and the expectation value of the current operator J = N L — N R = —i2d/da. 

To calculate the spin current we find the ground state wavefunction given by the eigen- 
state of the zero mode Hamiltonian. The zero mode Hamiltonian is obtained as a result of 
projecting out the non-zero mode field 9{x). This projection is obtained with the help of 
the Renormalization Group (R.G.) equations given in Ref. [14]. We consider the mesoscopic 
region, namely that the thermal length = hv s /ksT is larger than the circumference L 
of the ring. L defines the ring temperature T rmg = hv s /k B L. The backward coupling is 
renormalized and replaced by A e //. At a scale b we have: A (6) ~ Xb~ x where 1 < b < L/a 
and x = 2{2K — 1). For J± > J\\, x > and A(6) decreases. For J± ~ Jy we have x ~ and 
A(6) ~ A/(l + Xln(b)) decreases much more slowly. 



When J\\ > Jj_, X(b) grows with b at the scale b* ~ e 1 /^ ). When A(6) diverges, it signals 
that a spin gap is opened. Thus for a finite ring, gapless excitations are possible in the 
regime 7iJ± > J\\ > J± if the length of the ring L is shorter than exp(l/A), L < aexp(l/A), 
and the thermal length L T obeys L T > exp(l/A) > L. 

Once the particle hole boson 9{x) has been integrated away, one obtains a zero mode 
Hamiltonian. The zero mode eigenfunctions ip n>m (a, (5) of the Hamiltonian given in Eq. (16) 
have the form: 



^n,m(a,/3) = Ae t{5 - 1)a/2 e t{q - [q]W2 e ma/2 e tm ^ 2 



1 + C r e ira/2 

r=±2,±4,... 



(19) 



For simplicity we consider the case that q is an integer, q — [q] = 0. Under this condition 
the wavefunction in Eq. (19) is characterized by n = 0, ±1, ±2, ... and m = 0. 
The effective zero mode Hamiltonian for Q = has the form: 

Heff = ~2lf j2 + A e // cos ( 2a + 27r ^)' ( 20 ) 
where A e // = (J||/27r)F(L/a)[sin(27rg)/(27rg)]. Here F(L/a) is the scaling function 

F(L/a) = (L/a)~ x , x>0 and F(L/a) = (1 + Xln(L/a))-\ x ~ 0. 

Here, g measures the ratio of the Zeeman energy and the spin liquid lowest energy state 
excitation in the ring. When q <^ 1 (h — > 0) the last term in Eq. (20) is replaced by 
cos 2a. For large magnetic fields, q ~ 0.5 the cosine term can be ignored and the spin 
current is given by 1m — —9^b(v s /L)(1 — Q(9 )/'k). Using typical values of exchange energy 
J ~ 100K — lOOOK and magnetic fields of 5 — 10 Tesla one finds q ~ 1 justifying the neglect 
of the Ising term. In the remaining discussion we consider the modification of the current 
caused by the term A e //cos(2a + 2nq). 

Using Bloch theory we compute the eigenfunctions of the Hamiltonian in Eq. (20). We 
find ip n ,i( a ) = (l/V^r) exp(i(5— l)a/2) exp(ina/2)U n: i(a); n — 0, 1, 2, 3 are the Bloch states 
in the reduced zone and I is the energy band index. We restrict the discussion to the sector 
Q — 0. Therefore the possible states for J are n — 0, 2. Odd states are excluded since 
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they give states with Q ^ 0. Using Bloch states ^ n ,i{ a ) = ( a \ n ^) w ^h eigenvalues E n ^ we 
compute the spin current formula using Eq. (18) and find 

Im = -g^ EEi(«. Oe"^ 3, = -2^ B E E 7 e^-^. (21) 

^ n=0 J Z n=0 I Z \ an / 

Here Z is the zero mode partition function Z = X^ =0 J2i e~ En ' l ^ kBT and \n, I) are Bloch states 
with the eigenvalues E n j. The second part of Eq. (21) was obtained from Bloch's theory, 
see Ref. [15]. Using Eq. (21) we calculate the current for two cases: (a) weak coupling, 
\ e ff/(hirv s /2L) <C 1, and (b) strong coupling \ e ff/(hirv s /2L) 3> 1. 

(a) The weak coupling case: Because X e ff decreases as L~ x , x < 1 and the kinetic term 
(the term proportional to J 2 ) decreases as 1/L, the weak coupling region is achieved only 
for J ± 3> J|| or large magnetic fields | sin 2irq/2irq\ <C I. For this case we can use first 
order perturbation theory [16] to find the ground state Bloch wavefunction i^ n =o(a) for the 
Hamiltonian in Eq. (20): 

XeffL ( e~ i2a e i2nq e i2a e -i2nq > 



i(6-l)a/2 



1 + 



7tv s \2-(5-l) (5-1) + 2, 
Substituting Eq. (22) in the first part of Eq. (21) gives for the spin current 



(22) 



Im = -g^B-r 1 



L V 7T 



i _ 1 / J\\/J L \ 2 /L\ 2 ^ (sm2nq\ 2 ' 



(23) 



(4tt 3 ) 2 \1 - J\\/txJ±) \aJ V 27T( 1 J 
where v s = 2Jj_a(l — J\\/irJ ± ) and 1 > x > 0. Equation (23) is valid for values of J\\/J± 
satisfying (J\\/J±)\ sin2irq/2irq\ < (a/L) 1 ^ x . 

(b) The strong coupling case: For this case we compute the eigenfunctions of Eq. 
(20) using the "instanton" method given in Ref. [17] or the solution given on page 
231 in Ref. [15]. For large A e // the eigenvalues E n j = ± are given by E Ht ± = (hcu/2) ± 
(huj VD / ti) cos[27m + 2tt(5 — 1)], where the energy scale Tiu is determined by Tioj = 
(hvfy L)[(l/2n)(J\\/ J^L/a) 1 -^ sin2vrg/27rg|] with = J ± a/h and the band width 



is given by \f~D = [exp — uj J 27r y 2 cos(o; + 2nq)da\. Using the spin current given by Eq. (21) 
we find, 

Im = -^/i B ^sin[27r(l-fi(^o)/vr]v A D v /(8/7r)(Jj|/J ± )|s2n(27rg)/27rg|(L/a) 1 --. (24) 



From Eq. (24) we see that when J\\/J± increases the current is suppressed since the band 
width \f~D scales like exp(— yJ~J^fj±). The results in Eq. (23) and (24) are a consequence 
of the twist of the boundary condition obtained in Eq. (8c). In both cases the current 
was obtained for the even case, namely A N = 0. For the odd case A N = 1 and instead of 
1 — il(9 )/ir we have £l(9 )/n. 

Whether it is possible to detect the spin liquid current is an important question. Because 
of the Aharanov-Casher effect, when a metallic rod is inserted at the center of the ring one 
expects to find an electrostatic potential $e(0) ~ (1$ /L)(2ir) 2 , i.e., a voltage proportional 
to the spin current 1$ . For a magnetic field h ~ 10 Tesla, an exchange energy Jj_ ~ 2000 
Kelvin and a ring of 10 /xm, the Zeeman/exchange energy ratio is of the order of 10 allowing 
us to neglect the Ising term Jy and obtain a maximal current 1$ . For a ring of 10 //m 
we need to have a low temperature such that L T = TiVg/ksT > 10 /iin; corresponding to a 
temperature of one Kelvin. 

Materials that can be used to fabricate antiferromagnetic spin-1/2 ring structures include 
copper-benzoate. It is possible to fabricate such rings using plasma etching techniques in 
which a mesoscopic ring structure is etched out on a substrate and the desired magnetic 
material is deposited creating a magnetic ring structure. The persistent spin current can be 
measured using a scanning capacitance probe. Due to spin-orbit coupling the current will 
generate an electric field and therefore a potential difference. A small current in the ring 
will generate a voltage of the order of a nano-volt. An alternative approach is to attach 
two leads to a Luttinger liquid ring in the Aharonov-Bohm geometry with a crown shaped 
magnetic field. This generates a phase both in spin and charge. If the ring in the center 
is threaded with a flux tube [18] that exactly cancels the magnetic flux due to the crown 
shaped field, the remaining contribution is due to spin. Note that this flux tube will not 
affect the magnetic field on the ring. 

In conclusion, we have shown that a crown shaped magnetic field when acting on a spin- 
1/2 Heisenberg antiferromagnetic chain can give rise to a spin current proportional to the 
solid angle 0,(6 q) = 7r(l — cos#o)- The amplitude of the current is controlled by the ratio 
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J\\/J± and the Zeeman energy. A possible experimental realization has been suggested. 
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FIGURES 

FIG. 1. The geometry of a spin ring with circumference L in a crown shaped magnetic field 
h(x). The coordinate x is along the ring. The field has a tilt angle of 9q an d subtends a solid angle 
fi(0 o ). 
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